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Abstract. We consider three-dimensional models for rate- independent processes describing materials 
I undergoing phase transformations with heat transfer. The problem is formulated within the framework 

00 ' of generalized standard solids by the coupling of the momentum equilibrium equation and the flow rule 

I with the heat transfer equation. Under appropriate regularity assumptions on the initial data, we prove 

the existence a global solution for this thermodynamically consistent system, by using a fixed-point 
p I I argument combined with global energy estimates. 
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^ . 1 Introduction 
> 
00 

«) . Motivated by the study of shape memory alloys, we consider rate-independent processes describing 

' materials undergoing phase transformations. In the framework of generalized standard solids due to 

' Halphen and Nguyen (see |HaN75) ). the unknowns are the displacement field u and an internal vari- 

' able z and the problem is described by the momentum equilibrium equation combined with a flow 

, rule for the evolution of the internal variable. A very powerful tool to study such problems is the so 

called energetic formulation, introduced in |MiT041 IMieOS) and later on developed and intensively ap- 
plied in |FrM061 [MIRMI IMiR071 IMie071 iMiPOTl IMRS08] . Note that coupling rate-independent processes 
with rate-dependent processes makes, in general, the problems much more difficult; see for instance 
[EfMOGl [MPMQ81 IBaROSl IRouOQal iRouOQb) . 

In this paper, we are interested in coupling the rate-independent process with the thermal process. 



. which is not rate-independent, and viscous damping. The model is based on the Helmholtz free energy 

VF(e(u), z, Vz, 0), depending on the infinitesimal strain tensor e(u) == ^{'Vu+'Vu^) for the displacement 
u, where (•)^ denotes the transpose of the tensor, the internal variable z and the temperature 9. For 
simplicity, we will omit any dependence on the material point x E fl and t £ [0, T] with T > 0. We 
assume that W can be decomposed as follows 

W{e{u),z, Vz, 6) = Wi{e{u),z, Vz) - Wq{9) + eW2{e{u),z). (1.1) 

This partially linearized decomposition ensures that entropy separates the thermal and mechanical vari- 
ables (see (12.81) '). Let us emphasize that the last term 6W2{e{u), z) allows for coupling effects between 
the temperature and the internal variable, which is motivated by the phenomenological models for shape 
memory alloys presented in Section [2] Since coupling terms will appear both in the momentum equilib- 
rium equation and in the inclusion describing the evolution of the internal variable, this setting is more 
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general than the one presented in |RoulO) . We make the assumptions of small deformations. The problem 
is thus described by the following system 

- div(croi+Le(M)) ^ I, del = Dc(„) W^(e(u), z, Vz, 9), (1.2a) 
+Mi + cri„ 3 0, CTi„ = D^VK(e(u),z,Vz,e') -divDv^H/(e(u),z,Vz,6'), (1.2b) 
c{e)e - div(K(e(M), z, e)ye) ^ Le(u):e(M) + edtW2{e{u), z) + *(i) + Mz:z. (1.2c) 

Here ^' denotes the dissipation potential. As it is common in modeling hysteresis effect in mechanics, we 
assume that is positively homogeneous of degree 1, i.e., ^"(72) = 75'(z) for all 7 > 0. The viscosity 
tensors are denoted L and M, c{6) is the heat capacity and Hi[e{u),z,9) is the conductivity. As usual, 
('), and d denote the time derivative the i-th derivative with respect to z and the subdiffcrential 
in the sense of convex analysis (for more details see |Bre73| V respectively. Observe that (|1.2ap . (jl.2bp 
and (|1.2cp are usually called the momentum equilibrium equation, the flow rule and the heat-transfer 
equation, respectively. 

The paper is organized as follows. In Section [21 the thermodynamic consistency is justified and 
some illustrative examples are presented. Then the mathematical formulation of the problem in terms of 
displacement, internal variables and temperature is presented in Section [3l Our problem is reformulated 
in terms of enthalpy, which is a crucial ingredient to prove the existence result. Sections 21 [S] and [B] are 
devoted to the proof of a local existence result by a fixed point argument. More precisely, in Section |4l 
we consider first the system composed by the momentum equilibrium equation and flow rule for a given 
temperature 9 and we prove existence and regularity results. Next in Section |S1 we recall existence and 
regularity results for the enthalpy equation for any given right hand side. Then a local existence result 
follows in Section |B] by using a fixed-point argument. Finally a global energy estimate is established in 
Section [7] leading to a global existence result for the system (|1.2p . 



2 Mechanical model 

We justify here the thermodynamic consistency of the model (|1.2p . Starting from the Helmholtz free 
energy W ^ we introduce the specific entropy s via the Gibb's relation 

s = -I)eW{e{u),z,Vz,9), (2.1) 

and the internal energy 

Winie{u),z, \7z, 9) = W{e{u),z, Vz, 9) + Os. (2.2) 
Then the entropy equation is given by 

9s + div(j) = ^, (2.3) 
where j is the heat flux and ^ is the dissipation rate. We get 

^ = Le(u):e(u) -I- Mz:z + *(i) > 0, 

and, assuming Fourier's law for the temperature, we have 

j^-n{e{u),z,9)V9. 

We can check now that the second law of thermodynamics is satisfied if > 0. Indeed, assuming that 
the system is thermally isolated, we may divide (|2.3p by 9 and Green's formula yields 

Si\x^ [ div(«(e(«),z,9)Ve) | f hc{u):c{u)+Mz:i+'i'{z) 

Jn ^ Ja ^ 

^ 1 B(c(u),z,e)Ve-V9 ^^j^ I I.o(u):e(u)+Mz:z+'i!(z) ^ q 

Jfi Jn 

We differentiate now Win{e{u), z, Vz, 9) with respect to time, we obtain by using the chain rule and (|2.2p 
that 

W'in(e(u),z,Vz,0) =De(„)VK(e(M),z,Vz,0):e(it) 
+ 'D^W{e{u),z,Vz,9):z + 'D^^W{e{u),z,'^z,9)-Vz + 9s. 
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(2.5) 



We integrate (|2.4p over {I, thus we use the Green's formula and (|2.3p . we find 

Win(e{u),z,Vz,d)dx ^ / Dc(^u)W{e{u),z,Vz,d):e{u)dx + / D^W{e{u), z,Vz,9):zdx 
Jn Jn 

+ / D^^W{e{u),z,\7z,9)-Vzdx + / (div(K(e(u), z, 6')V0)+Lc(?i):e(?i)+Mi:i+*(i)) dx. 
Jn Jn 

On the one hand, we multiply (|f ■2ap by u, and we integrate this expression over fl to get 

Dc(u)W{e{u),z,Vz,9):e{u)dx + / Le(u):e(u)da; = / e-iidx. (2.6) 



On the other hand, the definition of the subdifferential 55' (i) leads to the variational equality associated 
to pr2b)) 

/ B^W{e{u),z,Vz,e):zdx+ / B^^W{e{u), z,Vz,9)-Vzdx + / Mz:zdx + / *(i)da; = 0. (2.7) 
Jn Jn Jn Jn 

We use (1^ and (lO) into (^3]) . we obtain 

Win{e{u),z,Vz,e)dx = / ^•uda;+ / K{e{u), z,9)\79-ridx. 
Jn Jon 

This means that the total energy balance can be expressed in terms of the internal energy, which is the 
sum of power of external load and heat. 
Note that from (|2.ip . we get 

s = DgWo{9)~W2{e{u),z), (2.8) 

and 

Win(e(u), z, Vz, 9) = Wiie{u), z, Vz, 9) + 9T>gW^{9) - Wo{9). (2.9) 

We use (|2.9I) into (I2.3p . we may deduce the heat-transfer equation (|1.2cp with the heat capacity given by 
c{9) = 9'DlW^{9). 

Motivated by the study of shape memory alloys, we will focus in the rest of the paper on the special 
case where Wi is given by 

Wi{e{u),z,Vz) = iE(e(u)-z):(e(u)-z) + f |Vz|2 + -ffi(z), 

where the internal variable z is a deviatoric d x d tensor. Hi is a hardening functional and the term 
■||Vz|^, 1/ > 0, takes into account some nonlocal interaction effect for the internal variable. Moreover we 
will assume 

VF2(e(u),z) = ati{eiu)) + H2iz). 

Here al, with a > and I the identity matrix, is the isothermal expansion tensor. Furthermore, the sum 
H{z,9) = Hi{z)+9H2{z) may be interpreted as a temperature-dependent hardening functional. Note 
that the system (|1.2I) is then rewritten as 

- div(E(e(u)-z)+Q6'I+Lc(M)) = £, 

S*(z) + Mi - E(e(M)-z) + B.Hiiz) + 9D^H2{z) - i/Az 9 0, 

c{9)9 - div(K(e(M), z, 9)V9) = Le(w):e(M) + 9{aiT:{e{u))+T>,H2{z):z) + ^{z) + Mz:z. 

We may illustrate the presented model by a nontrivial example, namely, a three-dimensional macroscopic 
phenomenological model for shape-memory polycrystalline materials undergoing phase transformations. 
This model was introduced by Souza et al ([SMZ98 ) and by Auricchio et al (see jAuP02j V The internal 
variable describes the inelastic part of the deformation due to the martensitic phase transformation and 
the hardening functional H takes the form 

i?(z,0)^='ci(0)Iz|+C2(0)|zp+x(z). 
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Here x ■ ^dcv ^ [0,+oo] denotes the indicator function of the baU {z e R"^^^ : \z\ < 03(6')} and the 
coefRcients Ci{6) are positive real numbers. Let us observe that ci{9) > is an activation threshold for 
initiation of martensitic phase transformations, C2{9) measures the occurrence of hardening with respect 
to the internal variable z and C3 (9) represents the maximum modulus of transformation strain that can 
be obtained by alignment of martensitic variants. 

The dependence of the coefficients Ci, « = 1, 2, 3, with respect to 9 is due to a strong thermo-mechanical 
constitutive coupling coming from the latent heat absorption or release, which is one of the main features 
of the behavior of shape-memory alloys (see |AuP04) for more details). 

For mathematical purposes, we should regularize the hardening functional as in |MiP07) . Namely, we 
replace H{z, 9) by 

H\z,9) c,{9)VS^+c,i9)\^\' + ^^%tt!n^\ 

where (5 > is a small parameter. Then, assuming that the mappings Ci, i — 1,2,3, are of class we 
can consider an afhne approximation of H^{z, 9) as Hf{z) + 6'-ff|(z) with 



H(iz) = c,^PT\W+C2\z\' + 



where Ci > 0, i — 1, 2, 3. 

Let us emphasize that the existence proof presented in the next sections can be easily extended to 
more general models of shape memory alloys for which Wi is given by 

Wi{eiu),z,Vz) = iE(e(u)-S(0)):(e(M)-^(z)) + f|Vz|2+Fi(z), 

where the internal variable z is a vector of M^^^, with N > 2, and E is an affine mapping from R^-i to 
the set oi d X d deviatoric tensors. In such models, N is the total number of phases, i.e., the austenite 
and all the variants of martensite, and the components zi, . . . , zat-i of z and zat == 1 — J2k=i are 
interpreted as phase fractions. Then E(z) is the effective transformation strain of the mixture, given by 

N-l JV-1 



E{z) ^fe^fc + (1- E ^''•)^^' 
fe=i fe=i 

where Ek is the transformation strain of the phase k and the temperature dependent hardening functional 
H{z,9) is the sum of a smooth part w{z,9) and the indicator function of the set [0, (see |MiT99[ 
IMieOOl IHaG021 IGMH02[ IMTL02[[GHH07] ). The hardening functional can be regularized in a similar way 
as in the previous example by replacing H{z,9) by 







k=l 



and we may consider an affine approximation of w{z,9) as wi{z) + 9w2{z). For more details on this 
example the reader is referred to |PaPllj . 



3 Mathematical formulation 

We consider a reference configuration £7 C K."^. We assume that is an bounded domain such that d^l 
is of class 0^+''. We will denote by R^ym the space of symmetric 3x3 tensors endowed with the natural 
scalar product v.w = tT(v~^w) and the corresponding norm \v\'^ = v:v for all v,w E Rsym- In particular, 
we assume that 

VFi(e(u),z, Vz) = iE(e(u)-z):(e(u)-z) + ^\Vz\^ + Hi{z) and W2{e{u),z) = ati{e{u)) + H2{z), 

where v > 0, a > 0, is the isotropic thermal expansion coefficient, E denotes the elastic tensor and Hi, 
i = 1,2, two hardening Junctionals. Given a function £ : ilx(0,T) — > M^, we look for a displacement u : 
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rix(0,r) ^ a matrix of internal variables 2; : Ox(0,T) -> R^^^ and a temperature : Qx{0,T) ^ R 
satisfying the following system: 

- div(E(e(M)-z)+a6'I+Le(M)) = £, (3.1a) 
d^{z) + Mi - E{e{u)-z) + B,Hi{z) + eV)^H2{z) - i^Az 3 0, (3.1b) 
c{9)0 - div(K(e(M), z, 9)\79) = Lc(u):c(u) + 6'(atr(e(u))+D^i72(z):i) + *(i) + Mz:z. (3.1c) 

We have naturally to prescribe initial conditions for the displacement, the internal variables, and the 
temperature, namely 

w(-,0) = u", z(-,0) = z", 9{-,0)^e°. (3.2) 
The problem is to be completed with boundary conditions. More precisely, we suppose here that 

where rj denotes the outward normal to the boundary dQ of The original problem (13. ip can be rewritten 
in terms of enthalpy instead of temperature by employing the so-called enthalpy transformation 

g{9) = d= / c{s)ds. (3.4) 
Jo 

Clearly, g is the unique primitive of the function c, which is supposed to be continuous, such that 17(0) ~ 0. 
Furthermore, we will assume that for all s > 0, c(s) > c*^ > where c'^ is a constant. Hence we deduce 
that g is a, bijection from [0, 00) into [0, 00). We define 



otherwise, 



«:^(e(z.),z,^)^=*-^^Mg^, (3.5b) 

where g~^ is the inverse of g. For more details on the enthalpy transformation, the reader is referred to 
|Rou09bj and the references therein. Therefore the system (13. ip is transformed into the following form 

- div(E(e(u)-z)+aC(t?)I+Le(M)) = £, (3.6a) 
a*(i) + Mi - E(e(u)-z) + D^ffi(z) + CWl),H2{z) - lyAz 3 0, (3.6b) 
d ~ div(K=(e(u), z, i9)V?9) = Le(M):e(u) + C{^){atr{e{u))+D,H2iz):z) + *(i) + Mi:i, (3.6c) 

with boimdary conditions 

^aa^O, Vz-7;|,„=0, K^Vi9-77|,„ = 0, (3.7) 

and initial conditions 

it(-,0) = u°, z(-,0) = z°, 7?(-,0) = i?° =5(61°). (3.8) 

The identity p.6cl) is called the enthalpy equation. As usual Korn's inequality will play a role in the 
mathematical analysis developed in the next sections. We have assumed that dfl is of class C'^'^^, so that 
the Korn's inequality holds, i.e. 

aC^^-" > V« e Rlin) : \\e{u)\\l^a) > C^^°™ll«llHMn)> (3-9) 

(see |KoU88llDTlL76] l. 

Let us introduce now the assumptions on the dissipation potential 5*, on the hardening functions Hi, 
i = 1,2, and on the data E, L, M, £, c c{9) and k = K(e(u), z, 9), which will allow us to obtain some 
regularity properties that are needed to prove the existence result. 

We assume that the dissipation potential 4* is positively homogeneous of degree 1 and satisfies the 
triangle inequality, namely, we have 

V7 > Vz e R3x3 . ^ l^{,z), (3.10a) 

3C* > Vz e Rl^^^ : < *(z) < C*|z|, (3.10b) 
Vzi, Z2 e Rl^,l : ^'(zi+Z2) < *(zi) + *(Z2). (3.10c) 
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It is clear that p.lOap . (jS.lObp and p.lOcp imply that is convex and continuous. We impose that the 
hardening functionals Hi, i — 1,2, belong to C^(R^g^^;M) and satisfy the following inequalities 



3c"\c"' > Vz £ Rl^^ : Hi{z) > c"^z\^-c"', (3.11a) 

3C^,' > Vz e R3x^3 . \j^lH,{z)\ < C^;. (3.11b) 

Note that (|3.11bD leads to 

3Cf > Vz G Rl^^ : |D,ff,(z)| < (l+|z|), \H,{z)\ < (l+|z|'). (3.12) 
The elastic tensor E : — s> £(Rg^j^, Rg^j,"^) is a symmetric positive definite operator such that 

3c^ > Vz e L2(l];R3yX^) : c'^||z||^2(o) < J Ez:zdx, (3.13a) 

V*,j,fc = 1,2,3: E(.),^M eL-(f]). (3.13b) 
We assume that L and M are symmetric positive definite tensors. This implies that 

3c^, > Vz e R^''^ : c^|zp < Lz:z < C^\z\^, (3.14a) 

3c^, > Vz G M3x3 . ^M|^|2 < j^^.^ < pM|^|2^ (3^^4,^) 

We assume furthermore that 

Vz e R3^^^3 . ^ ]r3x^3^ (3^^5) 

We consider that £ is an external loading satisfying 

e eIl\0,T;L'^{n)). (3.16) 

Finally, for the heat capacity c and the conductivity k"^ we assume that 

c : [0,oo) — > [0,cxd) is continuous, (3.17a) 

> 2 Bc^^ > > : c'=(l+6l)'3i-i < c(6l), (3.17b) 

k'^ : Mgy^xR^^^^^ xR Rgym is continuous, bounded and uniformly positive definite, i.e (3.17c) 

3c^^ > V(e,z,i?) e M-l^^ X R^^^^^J* x M V-y e R^ : K'=(e, z, ^9)wt; > c''>p, (3.17d) 

36"'" >0V(e,z,i9) GR^y"^^ xR3x^3 ^jg. [^'^(e, z, 7?)| < C"=\ (3.17e) 



Let us emphasize that the three dimensional model for shape memory alloys presented in Section [2] fulfills 
the previous assumptions so that we may apply the abstract result obtained in the next sections. 

Let us end this section by some comments about the proof strategy. In order to obtain a local 
existence result for the coupled problem p.Gp ^ p.Sp . a fixed point argument will be used. More pre- 
cisely, for any given we define 6 Ci"^) ^'^d we solve first the system composed by the momentum 
equilibrium equation and the flow rule p.lap - (j3.1bp . then we solve the enthalpy equation (I3.6cl) with 
^c^gj|y-j^ CW)- This allows us to define a mapping 

: i9 z9, 

and our aim is to prove that this mapping satisfies the assumptions of Schauder's fixed point theorem. 
Therefore, let us consider a given d G L«(0, T; LP(f])) with p > 1 and g > 1. We define 9 = C0)- Since 
C is a Lipschitz continuous mapping from R to R, we infer that the mapping 

01 : L«(0,r;LP(l])) ^ L«(0, T; LP(1])) 

is also Lipschitz continuous. Furthermore p.l7bp implies that 

ye G [0,oo) : |l((l+0)ft-l) ^=^gl(0) < g{0). 



Global existence result for phase transformations with beat transfer in shape memory alloys 



7 



Thus we have 

e [0,oo) : < cm < CiW = 5r'w, 

and ^ 

Vi? e R : |C(i9)| < + - 1, 

with d+ = max(7?, 0) for aU i9 e M. It follows that 

e Vi9 e M : \C{^)\ < - 1 < (§-1?+)^ (3.18) 

Hence, for all /3 e and for ah i9 e L9(0, T; LP(17)), we have = C(^?) £ L'39(0, T; L^P(rj)) with 

- ~ i 

In the rest of the paper, we will assume that q > A and p = 2. 

When there is not any confusion, we will use simply the notation X(f2) instead of X(f2; Y) where X is 
a functional space and Y is a vectorial space. 

4 Existence and regularity results for the system composed by 
the momentum equilibrium equation and the flow rule 

This section is devoted to the proof of existence and uniqueness results for the system composed by the 
momentum equilibrium equation and the flow rule p.lap - (j3.1bl) under the consideration that = C(i^) 
is given in a bounded subset of L''(0, T; LP{n)) with q — j3iq and p £ [4, min(/3ip, 6)] . More precisely, we 
look for a solution of the problem (P„z): 

- div(E(e(u)-z)+a6'I+Lc(?i)) = £, (4.1a) 
a^'(z) + Mi - E(e(u)-z) + B^Hiiz) + eV>^H2{z) - vAz 3 0, (4.1b) 

with initial conditions 

u{;0) = u°, z(-,0) = z°, (4.2) 

and boundary conditions 

"|an=0> V2-r,|,„=0. (4.3) 

Furthermore we will establish some a priori estimates and some regularity results for the solution of 

(Pn.). 

As a first step, we use classical results for Partial Differential Equations (PDE) and Ordinary Differ- 
ential Equations (ODE) to obtain an existence result. Let A : H^(ri) — > (H^(r2))' be the linear continuous 
mapping defined by 

V(u,«) e (Hi(f]))2 : (^M,v}(HMJ^))',HMo) - / M-^^u-.^vdz. 

Jn 

Classical results about elliptic operators implies that A generates an analytic semigroup on L^(f2), which 
extends to a C°-semigroup of contractions on LP(f2). We denote by Ap (resp. Ae) the realization of its 
generator in LP{fl) (resp. LP^'^{fl)) and by Xq p{Q) the intersection of interpolation spaces 

X,,p{n) ^^*-(LP(J7),2?(^p))i_|,, n(L^'/2(17),2?(^E))i_i,„ 

(see for instance |HiR081 iPrSOl] and the references therein). Here T>{Ap) (resp. 'D{Ae)) denotes the 
domain of (resp. Ap). 

In the sequel, the notations for the constants introduced in the proofs are valid only in the proof and 
we also use the set Qr x (0, r) with t e [0, T]. 

Theorem 4.1 (Existence for (Puz)) Let 9 be given in L9(0,r;LP(n)). Assume that u° € Hi(f]), 
z° G Hi(f^) and that ((3?T0)) . ((37T4)) and ((3?T6)) hold. Then the problem (|4J|) -([ilS ]) admits a 

solution lu,z) e Hi(0,r;Hj(17) x L2(17)) n L°°(0, T; Hj(f7) xHi(f7)). Furthermore if z° e Xg,p(f]) then 
z e L9(0,T;H2(17)) nC"([0,r],Hi(^^)) and z eLi{0,T;L^{n)). 
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Proof. Observe first that for all / £ L2(0, T; (Hi(r2); R^)') and for all u* e H1(17;R3) the following 
problem 

-div(Ee(M)+Le(u)) = /, 

with initial conditions 
and boundary conditions 

"Ian = 0, 

possesses a unique solution u = C{u* J) G H1(0, T; Hj(r2)) n C°([0, T]; Hj(17)). Note that 

V/i, /2 e L2(0, T; (H1(1]); R3)') Vu* e H1(17; M^) ; = C{u* , h) + £(0, h). 

and the mapping 

/:o(-) = /:(0, •) : L2(0, T; (HJ(1])')) ^ Hi(0, T; HJ(1])) n C0([0, T]; Hj(f])), 

/ '-^ w, 

is linear and continuous. Moreover, the classical energy estimate gives 

VtG[0,r]: <^\\e{u{-mlHu)+<^C''°''' (4.4) 

for aU / G L2(0,r; (Hj(17)')) and u ^ Co{f). It follows that and dlTTbl can be rewritten as 

+ Mi + + D^iJi(z) + eD,H2{z) - i^Az + 51(6*) + (72(^2) 3 0, (4.5) 

with initial conditions 

zi;0)^z"eil\n;Rl-J), (4.6) 

and boundary conditions 

= 0, (4.7) 
where 51(61) -Ee{C{u°,e)) - Ee(/:o(div(a6lI))) e Hi(0, T; L2(0; R^xS)) 

z ^ Ee(£o(div(Ez))). 

Let ip : L^{n;Rl^^) (-00, +00] defined by 



+00 otherwise, 



is a proper, convex lower semicontinuous function. Clearly, dtp is a maximal monotone operator, for more 
details on the maximal monotone operators and their properties, the reader is referred to |Bre73] . The 
resolvant of the subdifferential dip is defined by 

Ve>0: Je = {l+edipr\ 

Note that the resolvant J',: is a contraction defined on all L^(r2; R|^^^^). Let us introduce also 

Ve > Vz G L2(f];R3^^3) : ^,(z) min +<^(^)}, 

which, is a convex and Frechet differentiable mapping from L2(f];R3^^3) to R. Using |Bre731 Prop. 2.11], 
we know that the Yosida approximation of dip coincide with the Frechet differential of ip^, i.e., we have 

Ve>0: a^/='-i(I-J,), (4.8) 
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(4.11) 



and dipe is a monotone and i-Lipschitz continuous mapping on L^(f2; R^^^^). 
We approximate now the problem (I4.5I) - (|4.7I) by 

d^iQ + Mi, + dip.iz,) + Ez, + D,i7i(z,) 3 -{gi{e,)+g2iJeZe)+e,BzH2iJeZ,)), (4.9) 

with the initial condition 

z,(0,-) = ;2°. (4.10) 

Here 9^ £ Cg°(0, T)(g)Cg°(rJ) and d'i>{z^) is taken in the sense of the L^(ri)-extension of the subdifferential 
of the convex function (see examples 2.1.3 and 2.3.3 in IBre73] ). i.e., (|4.9p is equivalent to 

- g2{JMt,x)) - e,{t,x)V>,H2{JeZe{t,x))) £ d^{z,{t,x)). 

for all t G [0,r] and almost every a; S fi. Since we are looking for a solution with values in Kj^v^, the 
test-functions should be taken in Kj^^^'^. More precisely, we have 

«-(?) - *(ie(t,a;)) + (Mi^{t,x)+dip^{z^{t,x))^'Ez,{t,x)+T),Hi{z^{t,x)) 
+gi{e,{t, x))+g2iJeZe{t, x))+e,{t,x)D,H2{JMt, a;))):(z-z,(t, a;)) > 0, 

for all z G R^cv ^-'^'^ ^ ^ [0, 2^] and almost every x £ Q. Observe that (|4.1ip is equivalent to 

^'dov(z) - «'dcv(i£(i,a;)) +Projg3x3(Mi,(i,a;)+a<y9e(ze(t,a;))+Eze(i,a;)+D^i7i(ze(t,a;)) 

dcv 

+gi{e,{t, x))+g2{JeZ,){t, x)+0,{t, x)T),H2{JeZ,(t, x))):{z-z,(t, x)) > 0, 

for aU z £ R^^^ and t £ [0,T] and almost every x £ Q. Here ^dcv is the restriction of * to R^^^ and 
Projg3x3 denotes the projection on R^^^ relatively to the inner product of Mgy^. Using the definition of 
the subdifferential d'$dcv{Zi{t, x)) leads to the following differential inclusion 

- Projjg,3^x3{dLpeizeit, x))+Ez^{t, x)+'D^Hi{zt{t, x)) + gi{0t{t,x)) 

+g2iJMt, x)) + d,{t, x)^,H2{JeZe{t, x))) £ a*dov(i.(t, x)) + Mz,{t, x) , 

for all t £ [0,T] and almost every x £ ^. Note that the operator 9^dcv + M is strongly monotone on 
^dCT weU as its L2(rj)-extension, stih denoted 9^'dcv + M, on li^i^] H^dov )• This operator is invertible 
and its inverse is ^-Lipschitz continuous. Thus (|4.9p can be rewritten as 

i, = $(0„z,), (4.12) 

with 

$(0e,2) - (a*dcv+M)-i(ProjK3x3(-a(^,(z)-Ez-D,i7i(z)-gi(0,)-g2(Je;2)-eeD,i/2(Je2))), 

dcv 

for all z £ L^(0, T; L^(ri; Mdcv')) ^-"^^ we have to solve this differential equation for the unknown function 
z^. We may observe that we are not in the classical framework of ODE in L^(f7; Kd^ ) (^^^ |Car90] ). since 
g2{Jez) can not be defined for z £ \?{n-R^^^) but only for z £ h'^{{),T-\?{n-R^^^)). Nevertheless, 
we can solve (|4.12p with the Picard's iterations technique. Indeed, <^{9t,z) £ C°([0, T]; L^(ri)) for all 
z £ C°([0,T];L2(17)) and £ {[Q ,T]\\? [Vl)) is a solution of (gUl-gini) if and only if z, is a fixed 
point of the mapping 

A,: C'^([Q,T]-h\n))^G\[Q,T];h\n)) 

z^ K^{z) -A^ z^ + I $(6le(-,s),z(-,s))ds. 
Jo 

Let us assume that zi,Z2 £ C°([0, T]; L2(f2)). By using (I3.11bp and (|3.13bp . we find that 

||Ae(zi(-,t))-Ae(z2(-,t))||L2(0) < / ( II ^-^e (^1 (-,«)) "^^^e (^2 (', s) ) || (n) 

+ (||E||L^(O)+Cf/+Cf/||0e||L-(Q^))lkl(^s)~Z2(^s)||L^(O) + ||52(JeZl(-,s))-52(J.Z2(^s))|l^ 
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for all t G [0,r]. Since dipe is i-Lipschitz continuous on L^(f2) and 52 is linear, it follows that 
V< e [0,r] : ||A,(zi(-,i))~A,(z2(-,0)llL^(J^) < / ||zi(-,s)-Z2(-,s)||L2(n)ds 



(4.13) 

+ ^ / ll32(Je2^l(-,s)-J,Z2(-,s))||L2(n)ds, 



where C = ^ (i + ||E||Loo(n)+Cf^i +Cj^2 ||6'e l|L==(Qr)) • Thus the energy estimate (03) allows us to infer 
that 

Vs e [0,T] Vz e L2(0,T;L2(f7)) : \\g2{z{; s))h.^n) < j'^w^L ( / Viv(lE^(-, ^)) ||^h;(o))' dci) ' 



which, leads to 

V< G [0,T] : f \\g2iJM;s)-JM;smmn)ds 

< 4ife^ r{r\\jM:'y)-JM;^)\\hmd<yyds. 



Since J7e is a contraction on L^(il) we infer 
Vte[0,r]: ^'||g2(Je^i(^.)-J.^2(^.))||L^(o)d.< 

^ l|Ii-lll.°°(S"!) , II 11 

^ VcEc'-CKorn H I ^1 " ^2 1 1 C" (O , T ; (O) ) ■ 

Finally, we find 

Vt £ [0,T] : ||A,(zi(,t))-A,(z2(-,i))||L^(o) < j5fei£)i|Ni-^2||co(o,T;L^(0)). 
We iterate the previous computation to get 

W e [0,r] Vm > 1 : ||Ar(zi(-,t))-Ar(z2(-,i))||L^(n) < (C^+i51l^f£)™£|ki-.2||co(o,T;L^(o)). 



It follows that there exists mp > 1 such that 

\ ^M, /^E„Lr7Korn / mn! ' 



and A™'' possesses a unique fixed point in C°([0, T]; L^(51)), which is also the unique fixed point of Ag. We 
may conclude that there exists a unique solution e C^([0, T];L'^{D,; Kjcv )) to the problem (|O|) - (|4.10p . 
Let us choose now a sequence {0e)e>Q such that 

6,^6 in L«(0,T;LP(17)). 

Since p,q >2, it follows that 

-J> 6* in L2(0,r;L2(l])). 
We will establish that there exists a subsequence of {Zf )f-^o, stiU denoted by (ze)e>o, such that 

z,^z in Hi(0,r;L2(r!;R3^^^3)) ^^^j^^ 



Zf 



z in L°°(0,T;H1(1];M3x^3)) ^^^^ 



where z is a solution of the problem (|4.5p - (l4.7p . To do so, we notice that 

e [0,T] : z,{-,t) = (5*dov+M)-i(ProjR3x3(g,(t, J,z,(-,t))-9(^,(z,(-,t))-Ez,(-,i)-D,i/i(ze(-,i)))), 
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where g,{t,JeZei-,t)) -{g^{e,{;t))+g2{JM-,t))+S,{-,t)^MJM;i)))- Define 

Vt G [0,r] : w,{-,t) = g,{t,J,z,{-,t)) - d^e{ze{:,t)) -Ez,{-,t) ^D,Hi{z,(-,t)). 

We have e C"([0, T]; L2(1]; RSyX^)) and 

Vt e [0,r] : w,i-,t) + d^,{z,{;t)) + Ez,{;t) + D,Hi{z,i-,t)) ^ g,{t,J,z,{;t)), (4.14a) 
Vt e [0,T] : = (5^'dev+M)-i(ProjR3x3 i))). (4.14b) 



In order to obtain a priori estimates, we multiply (|4.14ap by ie, and we integrate this expression over Qr 
to get 

/ 'We.z^dxdt+ / d(p^{zg):z^dxdt+ / Ez^:z^dxdt+ / D^_ffi(Z(;):icdxdt = / gf:{t, J'^Zf:):Zf:dxdt. 

We observe that /g 9(/?e(ze):i:e da;di = (^e(2e(-, t)) — <i5e(2°), /g 'DzHi{z^):z^dxdt — j^Hi{z^{- ,T))dx — 
J^Hi{z°)dx. Moreover, recalling that z^ G C°([0, T]; L2(f]; R^xli^^^ ^^^^ 



/ ?iie:i£dxdt= / Projjj3x3 (wg):i£ dxdf > c*** / 
jQr JQr "'O 



l^e||L2(a)di- 



Since z° € H^(ri; Mj^^^'') = 2?(v?), we have (pe(z") < ^{z'^) (for technical details, the reader is referred to 
[Bre73j ) and 



dt + (pe(ze(-,r)) + / i?i(Ze(-,T))d2; + i / Ez, (• , t) (• , t) dx 



2 / ll-'e||L2(0) ' 

/o 



<¥>(z°)+ / i7i(z'')dx + i ( Ez°:z°dx + ^ f \\g,{t,JM; 
Jn Ja Jo 



L2(0) dt. 



Furthermore 



Ve > Vz G L2(r!;M3-^3) . ^ |||a^,(f)||2,^^^ + ^(j^^) ^ ^Hz-J^zH^.^^^ + ^(J,z) 

Thus J,z G CO([0,T];Hi(r!;R3x^3)) - ^ cO([0, T]; L2(r!; M^^^J*)) and 



¥ / P.llL^(o)dt+^lk.(-,r)-J,z,(-,T)||2 +^(J,z,(.,r))+ / Hi(z,(., T))dx 
+ i / Ez,(-,T):z,(-,T)dx<(p(zO)+ / iJi(z°)dx + i / Ez'>:z'' dx + ^ f ||.g,(i, J,z,(-, t)) || 



lL2(n) 



dt. 



By using p.llap . ((XT^ and (I3.13ap . we infer that 



ll^e|lL2(o)di+ilNe(-,r)-J,z,(-,T)||^,(^^ +Ci||z,(-,r)||2,^^^^ 



+ C^\\\/{JM;r))\\h(n)<Co + ^ / ||g.(i,J.z.(-,t))||^.(^,)dt, 



(4.15) 



where Co = (^(z") + ^ i/i(z") da: + B'-f^i + i /j-jEz":zOdx and Ci = min(c^i+^, f). Since J, is a 
contraction and J7e0 = 0, we have 



Vz G L^(f^;ML'v ) : II^.^IIl^(O) - WJeZ-JMLHn) < 



Then we find 



f / P.llL^(n)di + Ci||J,z,(-,T)||2j,(j,) <Co+ / ||5.(t,J.z.(-,t))||^(f2)dt. 
Jo Jo 



(4.16) 
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For the bound on ge{t, JeZe{-,t)), we use p.l2p to get 

/ llff.(i,Je^e)||^2(f,)dt <3 / ||.9l(0.)||^2(o)dt + 3 / \\g2{JeZe)\\l.^^^dt 

, (4.17) 

+ 6(Cf^)M \m\h^n)dt + 6{C^n' m'UeZel'dxdt. 

Jo JQ-r 

By using the definition of mappings gi and g2, the first two terms on the right hand side of (|4.17p can 
be estimated. More precisely, it is plain that there exists C2 > such that 



l.9i(^e)||^2(o)di+ / ||52(J.^.)||^2(o)dt 

/ r r r ^ ^''^'^ 

The last term on the right hand side of (|4.17l) is estimated by using Holder's inequality and the continuous 
embedding H^(r2) ^ L''(f2). Namely, it follows that there exists C3 > such that 



/ \ee\'\JeZe\^dxdt<Cs f ||^?,||2,(^JJ,z,|||,(^)dt. 

JQ-r Jo 

We insert (I4TT8)) and (|4TT9| into (|4TT6)) . we obtain 

^ £\\ie\\l^n)^lt + C,\\JM-,r)\g.^n)<Co + 

c.+6(cf2)2 r , r^!l!(^£!2!^^^^^ii7^ ii2 

i 2?t / l|f£|lL2(f2) ™ + / g?" IM£2e|lHi(0) 



(4.19) 



(4.20) 



Since 6',: is bounded in L*(0, T; LP(r2)), with q > 2 and p > 4, we may deduce from Gronwall's lemma 
that JeZe is bounded in L°°(0, T; Hi(ri)), ge{-,JeZe) is bounded in L2(0, T; L2(17)) and z^ is bounded in 
Hi(0,T;L2(n)) n L°°(0,T;L2(f])), independently of e > 0. Therefore ([HT^ implies that D,Mi{ze) is 
bounded in L°°(0, T; L^(il)), independently of e > 0. Furthermore the definition of the subdifferential 
d'^dcvize) enables us to infer from (|4.14b[) that 

- / Mz,:wPdxdt+ f |K||L2(o)dt< / (*«+Ze)-«'(ie))dxdt < f ^'«)dxdt, 

jQ-r JO JQr JQt 

where = Proj„,3x3 (uje), and using (jS.lObp . we get 

dcv 

IKI|L=(0,T;L^(n)) < C7*V^+C^P.||L^(0,T;L^(n))- (4-21) 

Observe that z^ is bounded in L^(0, T; L^(i7)), independently of e > 0, allows us to conclude that 
is bounded in L^(0, T; L^(r2)), independently of e > 0. On the other hand, we multiply (|4.14ap by 
dipe{ze) and we integrate this result over Q^. Recalling that dLpe{ze) takes its values in L^(ri; Mj^^^''), the 
Cauchy-Schwarz's inequality, p. 121) and (|4.2ip give 

I|5<^.(^.)I|l2(o,T;L^(o)) < (C^+C*)v^+(Cf^ + ll]E|| L°°(O))||^e||L2(0,T;L2(n)) 

+ C**||ie||L2(0,T;L2(f2)) + ||.9e(-i ^e^e) l|L2(0,T;L2(n)) • 

Thus d(pe{ze) is bounded in L^(0, T; L^(r2)), independently of e > and finally We is bounded in 
L2(0,r;L2(f7)), independently of e > 0. Furthermore, there exists z e Hi(0,r;L2(rj)) nL°°(0,T;L2(f])), 
z G L°°(0, T; H^(Sl)), w,v G L^(0, T; L^(ri)) such that it is possible to extract subsequences, still denoted 
by Ze, J^eZe, and dipe{ze) satisfying the following convergences 

Ze-^z in Hi(0,r;L2(f7)) weak, 
Ze^z in L°°(0,r;L2(f7)) weak*, 
JeZe^I in L°°(0,T;Hi(r2)) weak*, 
We-^w in L2(0,T;L2(f^)) weak, 
dip,{z,)^v in L2(0,T;L2(rj)) weak. 



Global existence result for phase transformations with beat transfer in shape memory alloys 



13 



Moreover, using (I4.15|) there exists C4 > such that 

Vr e [0,r] : ||z,(.,t)-J,z,(.,t)||2,(^) < ae{l+\\g,{; J,z,)\\l^,^^.^,^^y^) , (4.22) 

which, allows us to infer that z = z. Since J^^ is a contraction on L^(il), it follows also that, possibly 
extracting another subsequence, still denoted by J^z^, we have 

Vre[2,6): J,z,^z in C"([0, T]; L'-(fl)). 

We may deduce from (|4.22p that 

z,^z in C''(0,r;L2(f})). 

Our aim now consists to pass to the limit in (|4.14[) . To do so, we observe that the mapping £0 is linear 
and continuous, the mappings g\ and (72 are also continuous, which, gives 

g^{Q,)^g^{Q) in \? ,T -,1? {^)) , 
92{JeZe)^g2{z) iu (0, T; (17)) . 

Since DzHi, i = 1, 2, is Lipschitz continuous, it follows that 

DzHi{z^) +Ez^ ^DzHi{z) +Ez in C" {[0,T];L^{n)), 
Vr6[2,6): DzH2{J,z,) ^ DzH2{z) in C°([0, T]; L'-(17)). 

Hence we consider r such that | + i = ^ to get 

e,BMJeZe)^eDzH2{z) in L\0,T;L^{n)). 

We can pass to the limit in all terms of (|4.14al) . we obtain 

W + V + D,i/i(z) +Ez^ -{gi{9)+g2{z)+9DzH2{z)). (4.23) 

It remains to pass to the limit in (I4.14bp and to prove that z solves (|4.5p - (l4.7p . Using (|4.10p and the 
strong convergence of {zg)fyo to z in C'^([0, T]; L^(f7)), we infer that z(0, •) = Then we prove that 
v{t,-) g dip{z{t,-)) for almost every t € [0,T], i.e, v G dip[z) where dip is identified to its L^(0,T)- 
extension. To do so, we use the classical results for maximal monotone operators (see [Bre73] V More 
precisely, since dipi{z^) G dip{J^eZt), it is sufficient to prove that 



limsup / {dipt[z^)):{J^Zt)dxdt < / v.zdxdt, 



(4.24) 



which is an immediate consequence of the convergence results for the subsequences {dipeize))e>o and 
{<JeZe)e>o- Hencc, using the definitions of w and ip, we may conclude that z is a solution of (|4.5p - (|4.7p if 
we can prove that Proj^sxa (?ii(i, •)) — Mz{t, •) S 9^'dov(2(i, •)) fo'' almost every t £ [0,T]. So we may use 
the same trick as previously and we only need to check that 

limsup / Projg3x3(we):iedxdt < / Projjj3x3(w):idxdt. (4-25) 

e^O Jq^ lev Jq^ dov 

But takes its values in M.^^^ and thus z takes also its values in M.^^^ ■ It follows that (|4.25p is equivalent 
to 



limsup / 

e->-0 Jg. 



limsup / We'-Zedxdt < / w.zdxdt. 

'Qt Jqt 



We compute /^^ w^:z^ dx dt and w.z dx dt from (|4.14aP and (|4.23p . then the convergence results 
obtained above imply that (|4.25p holds if and only if 



liminf / dpe{ze)'.Zedxdt > / dp{z):zdxdt. 
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We observe that J^^ diff{ze):Zedxdt — ipe{Zi{T)) — ip^{z") > (p{J^Ze{T))~ipe{z") a.nd f^^ d(p{z):zdxdt = 
(p{z{T)) — ip{z'^)- Recalling that z° € T>{ip), we get \\m^^Qt^^{z^) = f{z^) and the lower semicontinuity 
of allows us to conclude. This proves the existence result. 
Finally we observe that 

-ProjR3xj!(w) = PY0]^3^^j{gi{e)+g2{z)+¥.z) + T>,Hi{z) + eD,H2{z) + v, 

and, using the definition of the mappings gi, 52 and if 

-Proj„3x3(u;) = Proj„3x3(-E(e(M)-z)) + DzHi{z) + eDzH2{z) - vl\z. 

dcv dcv 

So we may rewrite (|4.1bp as follows 

i - vMr^ I\z = m-^ r . (4.26) 

with = Projjj3x3(E(e(u)-;z)) - D^i/i(z) - B\i^H2{z) ~ ii and ^ = Projg3x3(w) - Mi e a^'dcv(i)- 
With ([XTU)) we infer that 

VV- e 9*dcv(i) : |V'I<C'* a.e. (a;, e f7 x (0, T). 

Furthermore, since z e L°°(0, T; Hi(0)), we infer with ((31^ that D^Hi{z) belongs to L°°(0, T; LP(0)) 
for i = 1,2. Then it follows that M"!/^ belongs to L9(0, T; L2(f7)). Since z° G X,^p(f7), we may deduce 
from the maximal regularity result for parabolic systems that z € L'^(0, T; H^(f7)) n C''([0, T]; H^(ri)) and 
i e L9(0, T; L2(fJ)). We refer to |Dor93[ iHiROSl IPrSOlj and the references therein for more details on the 
maximal regularity result for parabolic systems and its consequences. □ 

Let us observe that here neither + M nor dip + D^iJi + 6DzH2 + Projjj3x3o(E+(7i+g2) are 

linear and self-adjoint operators of L^(r2; Mj^^^'^), so we can not use the ideas proposed in jCoV90l to 
prove uniqueness. The uniqueness result proved below relies on the boundedness assumption (|3.11bp for 
the hardening functional Hi combined with Gronwall's lemma. More precisely let hi G C^(Ri^^y''; M) be 
defined by 

Vz e : hiiz) Hiiz) - C"- |z|2, (4.27) 

with a real number C^^ > 0. Assumption (jS.llbp implies that there exists C''^ > such that 

Vzi,Z2 eKdcv : \I^zhiiziy-D,hi{z2)\ < C'''\zi-Z2\. (4.28) 

Proposition 4.2 (Uniqueness for (Puz)) Assume that 9 is given in L«(0, T; LP(0)), (|3TT0)) . (|31T|) . 
(|3T4|) and (|3TT6l) hold and u" G z" G X,,p(f2). Then the problem (|4?T|) - (|4T3l) admits a unique 

solution. 



Proof. Let ^1 = (mi,zi) and ^2 == (^2,2^2) be two solutions of (|4.1ap - (l4.1bp satisfying (14.21) and (|4.3p . 
With the results of the previous theorem we already know that {ui,Zi) G C"([0, T]; Hg(r2) x H^(fi)) and 
Az, G L«(0, T; L^in)), i = 1, 2. Define 



j{t) = \ I E((e(ui)-zi)-(e(w2)-2:2)):((e(Mi)-zi)-(e(M2)-Z2))da: 

^ (4.29) 
+ f / V(zi-Z2):V(zi-Z2)da; + C^i / |zi-Z2pdx 



for aU t G [0,r]. Since 

7(i) i / E((e(Mi)-zi)-(e(M2)-Z2)):((e(ui)-zi)-(e(u2)-Z2))da; 
Jn 

A{zi-Z2y{zi-Z2)dx + C"^ I |zi-Z2pda;, 



the mapping 7 is continuous on [0,T] and its derivative in the sense of distributions belongs to L^(0,T). 
Thus 7 is absolutely continuous on [0,r] and with assumptions p.l2p and p.l3ap combined with Korn's 
inequality, we infer that there exists a real number k > such that 

lit) > K{\\uii;t)-U2i;t)\\l,^,,^ + ||zi(.,t)-Z2(-,t)||^i(o)) Vt G [0,T]. (4.30) 
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On the one hand, recaUing that i95'dov(0 is a monotone operator, the Green's formula and (|4.3p enable 
us to deduce from (j4.1b[) that 



Vz, £ L'in:R'^^^') : / {Mz,-Eie{u,)-z,)+B,Hl{z,)+e'D,H2{z^))■.iz,~z,)dx 
n 

V / Azj-(ii-Zi)da; < / (*dcv(zi)-*dcv(i:i)) dx. 



(4.31) 



On the other hand, we multiply (I4.1ap by iii — Ui, we integrate this expression over 51, and with the help 
of the Green's formula together with (|4.3p , we obtain 



yu, e nlin) : {E{e{u,)-z,)+a0l+'Le{u.,)):{e{u,) -e{u,))dx = e-{u,-u,)dx. (4.32) 
Jn Jn 

Thus we add (li3T|) and we get 

Zi-V.(e{ui)-Zi)+'DzHi{zi)+9T):,H2{zi)):{zi-Zi)Ax - V / Azi-(ii-Zi) da: 
o Jn (433^ 

(E(e(Mi)-2:i)+aOT+Le(ui)):(e(wj)-e(ui))da; < / (^•(wi-Mi)+^'dev(2i)-*dev(ii)) dx. 



We choose now {ui,Zi) = (u^-i, z^-i) for i = 1,2 in (|4.33p . we add these two inequalities, and with the 
help of ([XTi]) . we find 

7(t)+cM / |ii-i2pdx + c^ / |e(iti)-e(w2)pda; 

Jn 

< - / (D,/^i(zi)-D,/ii(z2)):(ii-i2)da:- / 0(D,i/2(2i)-D,i/2(22)):(ii-i2) dx a.e. te[0,T]. 
Jn Jn 

We estimate the terms of the right hand side of the previous inequality as follows: 

/ {T),hi{zi)-Ty,hi{z2))-{zi~Z2)dx < C^^ / \zi-Z2\\zi-Z2\dx 
Jn Jn 

<f /|ii-i2pdx + i^ /|zi-Z2pda:, 
Jn Jn 

and 

/ 0(D,ff2(zi)-D,i?2(22)):(ii-i2)dx < Cf^^ / \e\\zi~Z2\\zi~Z2\dx 
Jn Jn 

<^ [ \i,-i,\^dx+^ [ \e\'\z,-Z2\'dx<^ [ \i,-Z2\'dx+^\\9\\l^^^\\z,-Z2\\l^^y 

Jn Jn Jn 

Since H^(r2) ^ L^(r2) with continuous embedding, we infer that there exists C > such that 

7(t) < C{l+\\e{-MlHn))\\^^i-^tyM-MHHn) < ^{^+Pi;t)\\lHn)Ht) a.e. t e [0,r]. 
But 9 e L9(0, T; LJ'(f7)) with g > 4 and p > 4, thus we get 

G [0,T] : 7(t) < 7(0)e^^«(^+"'(•■^)"^^<-))'^^ 
which, allows us to conclude. □ 

We provide that I? {u, z) is a continuous mapping from L«(0, T; LP(rj)) into Hi(0, T; Hj(f7) xL2(rj))n 
L°°(0,T;Hj(rj)xHi(ri)) where {u, z) is the unique solution of (Puz) when 9 = C(i^)- 

Lemma 4.3 Assume that ([3J0)) . ((3?TT|) . ((3ll| and ((3?T6| /loZd and f/iat m° G Hi(fi), z° e Xg.p(ll). Then 
-d^ {u,z) IS continuous from L«(0, T; LP(rj)) into H^O, T; Hj(17) xL2(ll))nL°°(0, T; Hj(17) xHi(17)). 
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Proof. We consider e L'?(0, T; LP{n)) and for i = 1,2, we denote by 6*, = C{^^) e L*(0, T; LP{n)) and 
(wi, Zi) the solutfon of the fohowing system: 

- div(E(e(u,)-z,;)+ae'J+Le(u,)) = £, (4.34a) 
d^{i,) + Mi, - E(c(u,)-z^) + D^iJi(z,) + 9B,H2iz^) - uAz, 3 0, (4.34b) 

together with initial conditions 

u,(0,-) = u°, z,(0,-) = ^°, (4.35) 

and boundary conditions 

Ui\3u = 0' ^Zt-V\oa = 0- (4-36) 
Since the mapping (f)i : -d ^ 6 ^ is continuous from L«(0, T; LP(f])) to L9(0, T; LP(f])), it is enough 
to check that the mapping 9 = C{d) {u, z) is continuous from L«(0, T; LP(17)) to L«(0, T; LP(f))). 

Once again the key tool is the structural decomposition (|4.27p combined with Gronwall's lemma. 
We reproduce the same kind of computations as in Proposition 14.21 More precisely, the total energy 
inequality associated to (|4.34p - (|4.36p is obtained by multiplying (|4.34ap and (|4.34bp by ii^-i—Ui and 
z^^i—Zi respectively, and integrating over 51. Then, we add these two inequalities, we find 



E(e(ui)-Zi):((e(u3_j)-i3_,)-(e(ui)-ii))da; + / Mi;j:(i3_.i-ii) dec 

Jn 

+ / a6'itr(e(u3_i)-e(wi))da:; + / Le(ui):(e(u3_,)-e(ui))da; 



(4.37) 

i'/\zi-{z^^i~Zi)Ax + / 'DzHi{zi):{z^^i~Zi)Ax 



+ / 6,'D^H2{zi):(zy,^i-z,)dx- / £-(u3-,-iiO dx + *dcv(i3-0 - *dcv(i^) > a.e. te[0,r]. 
Jn Jn 

It is convenient to introduce the notations: u = ui~U2, z = z\—zi and e = 91-92- Therefore, we take 
z = 1, 2 in (I4.37P and we add these two inequalities, we get 



E(e(?i)— z):(e(?i)— z) da; + / Mz:zd2;+ / Le(u):e(u) dx — / vAz-zdx 
! Jn Jn Jn 

+ [ {D,Hi{zi)-D,Hi{z2)):zdx < - [ a9tr{e{u))dx ~ [ {9i^,H2[zi)-92D,H2{z2)):idx. 
Jn Jn Jn 

Define 

V< e [0,T] : 7(<) = i / ¥.{{e{ui)-zi)~{e{u2)-Z2)):{{e{ui)-z{)-{e{u2)-Z2))Ax 
Jn 

-f / V(zi-Z2):V(zi-Z2)dx + C7^i / |zi-Z2pdx. 
Jn Jn 

As in Proposition I4.2[ we can check that 7 is absolutely continuous on [0,T] and (13.141) and (|4.27p imply 

that 

m+c'^\\-4hin)+c'^\\<^)\\lHn) < - j^{D.hi{zi)-B,hi{z2)):tAx 

- [ a9tT{e{tL))dx~ [ {9iT),H2{zi)-92DzH2{z2)):zdx, 
Jn Jn 

for almost every i e [0,r]. Clearly, it follows from (|4.28l) that 

7W+c^II^IIl^(o) +c^l|e(^)ll£.(n) <- J^^a9tTieiU))dx 

- [ {9iT>Mzi)~92V),H2{z2)):zdx + &'' ( \z\\t\dx. 
Jn Jn 

We estimate the first and third term on the right hand side with the help of Cauchy-Schwarz's inequality, 
while for the second term, we use the following decomposition 

{9iD,H2{zi)-92D,H2{z2)):i^{9D,H2{zi)+92{D,H2{zi)-D,H2{z2))):t. 
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Hence, we obtain 

at) + '-^mUn) + ^Mmhm < ^mUn) + &o\\hin) 

r (4-38) 

+ / {\e\\D,H2{zi)\\i\ + \e2\\B,H2{zi)-D,H2{z2)\\t\)dx. 
Jn 

It remains to estimate the last term on the right hand side of (|4.38p . We use (jS.llbp and (|3.12l) to get 

/ {\e\\B,H2izi)\\z\ + \02\\B,H2izi)-B,H2{z2)\\z\)dx 

<C^^ [ (l+|zi|)|^||l|da; + Cf^^ / \e2\\z\\t\dx. 
Jn Jn 

The Yomig's inequality imphes that there exists 7^ > 0, i = 1, 2, 3, such that 

^^(|0||D,ff2(2i)||^| + |e2||D.i72(^i)-D,iJ2(^2)||^|)da;<^p||2,(j^^ 

+^ / i^Tkpd.+^ / i^2n.td.+ ^"-<^-n^H^-^-^- iizii^.(,). 

Jn ' Jn 



Since zi G L''(0, T; H^(ri)) and ^ L°°{fl) with continuous embedding, we get 



Cf2 „;j„2 

L2(n) 



^ {\e\\B,H2{zi)\\t\ + \&2\\BMzi)-'D,H2{z2)\\z\)dx < g 



(4.39) 



We insert (|4.39p in (|4.38p and we choose 71 = 72 = j^pr^ ^nd 73 — ^^prj- Therefore the continuous 
embedding lV-{n) ^ L^{n) and KM give 

at) + ^mhin) + TM^)\\hin) < Cie,z,) + C(02)IN-||?ii(o) < C{e,z,) + ^7W, (4.40) 

dot .(3a)^_^ 2(Cf2)2 . - ^2(Cf2)2 - 



for almost every t G [0,T]. Here C(0>i) + ||0-||2^^^^ + iii^||^,||2^^^^||0-|j2^^^^ and 

C(^2) = — I H^2|Il4(o) where C > is the generic constant involved in the continuous 

embedding Hi(f7) ^ L'*(0). Since g > 3 we can check that C{e,zi) e h^iO,T) and 



C{d,Zi)dt< (^l^+^^'^Ja '' ) ll^llL2(0,T;L2(n)) + ^'■*^c« ^ ll^llP^, , ,Jl^llL«(0,T;L2(n))- 

L <!-2 (0jT;L°°(O)) 

Note that (7(6*2) G L^(0,r), which, thanks to Gronwall's lemma, gives 

Vte[0,r]: 7(t)< Ci0i;s),zi{;s))eifsC{ei{;a)-s(-,.))d.^^^ 



Recalling that the mapping ■& 1-^ 9 is Lipschitz continuous from L^(0, T; L^(r2)) to L''(0, T; L^(r2)) and 
maps any bounded subset of L'(0, T; L^(f2)) into a bounded subset of L''(0, T; L?'(r2)) with p = 2 and 
q > 4, p > 4, q = /Siq > 8, the last estimate allows us to conclude. □ 

Let us observe furthermore that the image of a bounded set of L'^(0, T; L^(i7)) by the mapping i? 1— >■ 
(0) (u,z) is a bounded subset of H^O, T; Hi^(f7) xL2(17)) n L°°(0, T; Hj(17) xH^ (rj)). 

Let us introduce now some new notations. For any r > 1, let 

V(r!;R3) = {u e L2(ri;R3) : Vu G L'"(fl; R3x3)| ^ 

and for any r > 2, let 
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We endowed V''(ri;R"^) with the foUowing norm 

The aim of the next two lemmas is to prove further regularity results for the solutions (u, z) of the 
system composed by the momentum equilibrium equation and the flow rule. More precisely, assuming that 
6 remains in a bounded subset of L''(0, T; LP{n)), we will prove that e{u), z and z remain in a bounded 
subset of L9(0, T; LP(f2)), L«(0, T; LP/^{n)) n L9/2(0, T; LP{n)) and L«(0, T; H2(rj)), respectively. 

Lemma 4.4 Assume that (|3.10p . (|3.11l) . (I3.13bp . (|3.14l) and (|3.16p /lo/d. ^sswme moreover that z° G 
X,.p(rj) and u° G Vg(fi;M3). Then e{u) belongs to W^'i{0,T]LP{n)) and ^ e(u) maps any bounded 
subset o/L9(0,T;LP(f7)) into a bounded subset o/ Wi'9(0,T;LP(r2)). 



Proof. The idea of the proof is to interpret (|4.1ap as an ODE for u in an appropriate Banach space. More 
precisely, let Jp = L^(ri; R^) xLP(r2; R^y„^) endowed with the norm 

It follows that T,p is a Banach space. Let us introduce now the mapping Ae 

Ae : Vg(r!;R3) ^ Jp, 

ip""^ (0,Ee(u)). 

Since E e L°°(ri), we infer that is a linear continuous mapping from VQ(r2) to J'p. Besides since L 
is a symmetric, positive definite tensor, classical results about PDE in Banach spaces imply that, for all 
(f = {(pi, (P2) S J^p, there exists a unique u € Voi^', R'^), denoted u = Ap{(p), such that 

yveV{n;R^): / Le(-u):e(i;) dx = / tpi-vdx+ / ip2:e{v)dx, 
Jn Jn Jn 

where p* is the conjugate of p, i.e ^ + i = 1. Furthermore there exists a real number C > 0, independent 
of ip, such that 

IMvpiam < C(ll</'i||L2(n) + ||</'2||Lp(o)) = ^11(^11^^, 

and Ap is linear continuous from Tp to Vq(17;R'^) (for more details, the reader is referred to |Val88] ) . It 
follows that (j4.1ap can be rewritten as 

u ^ Gp{iPze,u), (4-41) 

with ifzff = {£,Ez-a91) and Gpiipze, u) = Ap^ip^^g-AEu) = Ap{ipy^g)-Ap{AEu). 

With the assumption p.l6p and the previous results, we already know that (fze G L''(0,T;7^p) and 
we can apply classical results for ODE in Banach spaces to conclude that u S W^'''{0,T;Vq{Q;R^)), the 
reader is referred to |Car90j for more details. 

We can also obtain estimates for u and ii in VQ(il;R'^). To this aim, we introduce the following 
notations: Ca^ = \\K\\c{J^p,vl{n)) and C^j. = \\AE\\c{vf,{n),j^p)- Then, we observe that (I4.4ip gives 

M-,t)\\yp,^n)<CAp{\\^z0{-,t)~AEU°yp+CAj\u{-,t)-u"\\vp(n)) a.e. t e [0,T]. (4.42) 
Let us turn now to the term ||u(-, t)— M'^||vp(f2)- It is clear that 



\\u{-,t)-u llvp(n) < / l|5p(<^ze(-,s),u(-,s))||vp(n)ds 
Jo 

<Cap / \\iPze{-,s)-AEU°\\jr^ds + Ca^Cae / s)-u°||vp(f2) < 



Therefore, we may infer from Gronwall's lemma and the continuous embedding Il^(ri) ^ ljP{fl) that 
there exists a generic constant Ci > such that 

||^/(-,t)-MO||vP(n) <e^^pC'^^*CA^ / \\vze{syAEU°y,ds 

^ (4.43) 
<e^^p^^E*CA^Ci / (|K(s)||L2(o) + l|E||Lo=(o)ll^(-,s)llHi(o)+a||0(-,s)||Lp(n)+C'^.l|u°||vno))ds. 
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We insert (I4.43P in (|4.42l) . we find that there exists C*2 > such that 

ll'"l|L'!(0,T;V!'(n)) < <^2(|K||L<!(0,T;L2(n)) + ||IE||L~(n)|k||L<J(0,T;Hi(n))+a||^'||L9(0,T;LP(O)) + ||M°||vp(n))- 

Recahing Lemma [4.31 we infer that the image of any bounded set of L'?(0, T; LP(ri)) by the mappings 
d ^^ C{d) = 9 ^ e{u) and d ^ (0) ^ 9 ^ e(u) are still bounded sets in L°°(0, T; LP(f])) and 
L9(0,T;LP(r2)), respectively. □ 

Let us conclude this section with some regularity results for z and z. 

Lemma 4.5 Assume that ((3l0)) . ((3lT1) . (|3.13bp . ((3Jl) and ((3lB hold. A ssume moreover that z £ 
X,.p(17) and vP e Vg(f7;M3). Then z and C^z belong to L9/2(0, T; LP(r2))nL9(0, T; Lp/2(0)) and z e 
Co'([0,T],X,,p(fi))nL«(0,T;H2(rj)). Moreover B ^ (i,Az,z) maps any bounded subset of Li{0,T; LP (Q)) 
into a bounded subset of (Li/'^{0,T;LP{n))nLi{0,T;LP/^{Q))f x (C°([0, T]; X,,p(f7))nL9(0, T; H2(f7))). 

Proof. Let us rewrite again (j4.1bp as 

i - i/M^^Az = M"\r- (4.44) 

with = Projjj3x3(E(e(u)-z))-D^i7i(z)-6ID^i72(z)-^ and^ = Projjj3x3(w)-Mi e d^dcAz) where 
w has been defined in the proof of Theorem 14. II With the assumption (IS.lObp . we have 

VV' e ai'dcv(i) : |V'I<C'* a.e. {x,t) e Q x {0,T). 

Then Lemma 14.41 enables us to infer that Proj„3x3 (Ee(u)) — tp remains bounded in L°°(0, T; LP(ri)). 

dcv 

Furthermore, we know with Lemma 14.31 that z is bounded in L°°(0, T; H^(il)), so, using p.l2p . we 
infer that B^H^{z) is bounded in L°°(0, T; LP(f^)) for i = 1,2. Then it follows that is bounded 

in L'J(0,r;LP/2(f])) if g belongs to a bounded subset of L«(0, T; LP(f2)). We may deduce from the 
maximal regularity result for parabolic systems that z is bounded in L''(0, T; H^(ri)) and i is bounded in 
L9(0,T;LP/2(n)) (see |Dor93l[HiR08l IPrSOlj l. 

Since H^(ri) ^ L°°(r2) with continuous embedding, z is bounded in L''(0, T; L°°(f7)) and thus 
9T)^H2{z) is bounded in l.il'^{Q,T;lJ'{n)). We may deduce that is bounded in L«/2(0, T; LP(r2)) and, 
the maximal regularity result for parabolic systems allows us to infer that i and Az belong to a bounded 
subset of L9(0,r;LP/2(f7)) n L9/2(0, T; LP(r2)) and z belongs to a bounded subset of C0([0, T]; Xq,p(f7)) 
whenever 6 belongs to a bounded subset of L''(0, T; LP(il)). □ 



5 Existence and regularity results for the enthalpy equation 

In this section existence and uniqueness results for the enthalpy equation are recalled and some regularity 
results are obtained. More precisely, let us consider the enthalpy equation (P^): 

li* - div(K"Vi?) = / (5.1) 

together with initial conditions 

?9(0) = ?9°, (5.2) 

and boundary conditions 

Ti'V^-ri\,,, - 0. (5.3) 

We assume that the initial enthalpy belongs to L^(r2) and / belongs to L^(0, T\ L^(ri)). Furthermore, 
we assume that e L°°(QT;Ksym) and satisfies 

3c'^^ >0'iv ieW^ : K''{x,t)vv>c'^''\v\^ a.e. (x,t) Qt, (5.4a) 

3(7"' > : \ii''{x,t)\ < a.e. {x,t) G Qt- (5.4b) 
The weak formulation of the problem is given by 

Find d : [0,T] ^ such that ^ and for ah ^ G H^^^), 

i^^dx^- / 'K'^V'd-'\/£,dx — i fidx in the sense of distributions. 



(5.5) 
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Theorem 5.1 (Existence and uniqueness for (P^)) Under the previous assumptions, the problem 
(|ET|) - (iO)) possesses a unique solution d G C°{[0,T];'L^{n))nL^{0,T;ll^{n)) with d £ L^{0,T; (Hi(17)'). 
Moreover we have 

VtG [0,T] : ||^(r)||2^(,^,+2c«^ ri|Vi?W||£.(0) dt < e-(||i?°||^.(^) + ||/||2^(Q ). 

JO 

Proof. The proof of existence and uniqueness of a solution is quite classical and can be found in [Bre83i 
IRaT83j . The estimate is straightforward and its verification is left to the reader. □ 

Let us introduce the following functional space 

W = {d G L^{0,T;ll\n))nL°°{0,T;L^{n)) : ^ G L2(0,r; (Hi(f7)')}, 
endowed with the norm 

Vi9 G W : \\d\\w = lli?IlL2(o,T;Hi(n)) + \W\\L°-(o.T;L^n}) + \\'^\\h^o,T:{H^ny)- 

Due to |Sim87| . we know that W is compactly embedded in L''(0, T; L^(ri)). Note that the previous 
estimate implies that there exists a generic constant C > such that the solution of problem (P^) 
satisfies 

\\nw<C{\\i9'^h2^n) + \\fhHo,T;mn)))- 



6 Local existence result 

This section is dedicated to the proof of a local existence result for p.6p - p.8p by using a fixed-point 
argument. To this aim, for any given ■& G L^(0, T; LP(17)), we consider k"^ = k"^ {e{u) , z , 0) and / = 
= Le(u):e(u) + 6'(atr(e(M)) + D^H2{z):z) + + Uz:z in (P^), where (u, z) are the solutions of 
(Puz) with 9 — ({i}). With the results obtained in the Section lU we already know that belongs to 
L'?/4(o, T; LP^^{n)). Since p > 4 and g > 8, we infer that belongs to L^{0, T; L^{n)) and we can define 
"0 G C°([0, T]; L^(i7)) n W as the unique solution of (Pij). Thus we can introduce the fixed point mapping 
(t):^^T} from L'?(0, T; LP{n)) to L^O, T; LP{n)). 

Proposition 6.1 The mapping (p is continuous from L«(0, T; L?^(ri)) to L'?(0, T; LP(ri)). 

Proof. Let (i?„)„gN be a converging sequence of L^(0, T; L^(r2)) and let i9* be its limit. We denote by 
^„ 0(i?„) for ah n > and = Since (i?„)„gN is a bounded family of L9(0, T; LP(r2)), we 

infer from the previous results that ('i^„)„gN is bounded in C°([0, T]; L^(r2)) n W. We may deduce that 
(i?„)„gN is relatively compact in L^(0, T; LP(r2)) (see |Sim87j ) and there exists a subsequence, still denoted 
{'&n)neN, such that 

i?„^i9 in L2(0,T;H^(rj)) weak, 
i?„->i9 in L«(0,T;LP(17)). 

Let us define Vt = {w G C°°([0,T]) : w{T) = O}. Hence we observe that for all n > 0, we have 
G Hi(17) Vw G Vt : - /" dn{x,t)^{x)w{t)dxdt+ f K';ydn{x,t)V£.{x)w{t)dxdt 

Jqt ^ JQt (g^^^ 

f'^"{x,t)^{x)w{t)dxdt+ / ^°{x)^{x)w{0)dx, 



Qt -io. 

with kJ^ K'^(e(tt„), z„, 6'„) and (u„,z„) solutions of (Puz) with 9n Ci'&n)- Since (i?ri)nGN converges to 
•di, in L'(0, T; L^(r2)), we infer from Lemma l4.3l that (u„, 2;„)„gN converges to the solution (u*, z*) of (Puz) 
with J** = C0*)J-o^ Hi(0,T;Hi(17)xL2(17)) n L°°(0,T;Hi(rj)xHi(rj)). Let us recah that the mapping 
(j)i:-d^e^ C{-d) is Lipschitz continuous from L«(0, T; LP(r2)) to L«(0, T; LP(rj)), which implies, possibly 
after extracting another subsequence, that 

0n,u„,Zn ^ 0*,u^,z^ a.c. {x,t)eQT. 
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Note that the continuity of the mapping k"^ gives 

= K''(e(M„), z„, 0„) ^ = K''(e(M*), z*, 0*) a.e. {x,t) £ Qt, 
and due to the boundedness assumption on k"^, we obtain with the Lebesgue's theorem that 

K^V^w-^K^V^w in L^{0,T;L^{n)). 
Therefore it is possible to pass to the hmit in ah the terms of the left hand side of (|6.1|) to get 

e tf(rj) Vu; G Vt : - / i){x,t)^{x)w{t)dxdt + [ K';V'd{x,t)V£.{x)w{t)dxdt 
Jqt JQt 

= lim / f--{x,t)£^{x)w{t)dxdt+ [ i}°{x)^{x)w{0)dx. 

Recalling that p € [4, 6] and <? > 8, we infer that the mapping n> is continuous from L^(0, T; L^(i7)) 
to U^{0,T;U^{n)) with 77 = f + 5 and ^ = i + i. Indeed, for any in L«(0, T; LP(f7)), let (w„ z^) be 
the solution of (Puz) with 9i = ({"di), i — 1,2, we find 

_ f2 ^ Le(ui+W2):c(ui-it2) + (0l-^^2)(atr(e(ul))+D,i^2(zl):il) 
+ 92{atT{c{ui-U2))+'DzH2{zi):zi-I),H2{z2)-Z2) + ^'(ii) - ^'(^2) + M(ii+i2):(ii-i2)- 

On the other hand, p.lOcP and (|3.10bp imply 

|*(il)-*(i2)| <C*|ii-i2|, 

and ([XT2|) and (|3.11b|) give 

\D,H2{ziy.Zi-D,H2{z2):Z2\ < \B,H2{zi)\\zi-Z2\ + \D,H2{zi)-D,H2{z2)\\z2\ 

< C^'il + \Z1\)\ZI-Z2\ + C^z'\zi-Z2\\z2\. 

The boundedness and the continuity properties proved in Lemma and Lemma respectively, allow 
us to deduce the desired result. Therefore, we may infer that 

£ R\n) \fw e Vt ■■ lim / f^^{x,t)^{x)w{t)dxdt=f f^' {x,t)^ix)wit)dxdt. 

We conclude that i? is solution of problem (P^) with the data and /''* . Moreover by uniqueness of 
the solution, it follows that ^ — -d^, and the whole sequence {'&n)nefi converges to i?* = C(^*)- O 

We establish now that the mapping (p fulfills the other assumptions of the Schauder's fixed point 
theorem. To this aim, we introduce some notations: let R'^,R^ > be any given positive real numbers 
such that max(||M°||vP(n), lk°||x,,p(n)) < R° and ||i^||L4(o,T;Lp(n)) < R"^ ■ Clearly, we have 

IICWIIl,(0.T;L.(O)) = ||0||l,(o,T;L.(o)) < {^)^m'^\mI}io,T-,LPm<R' = i^^^^ 

Therefore once again the results of Section |4] are used, which imply that there exists a constant R^ = 
Rf{R°,R^, |K||c"([o,T]:L2(o))) > 0, depending only on R^ and \\i\\co{[Q,T]■,L^n)), such that 

ll/''llL<!/4(0,T;LP/2(n)) < R-'^iR^^R^, ||^||c0([0,T];L2(O)))- 

The results of Section [S] imply that there exists a generic constant C > such that 

< c(|^^|^^'^i^^(i^°,i^^ m\coi[o,nL-m))+\\^''hHn)). 
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Now let < T < T and let us introduce the following functional space 

Wr = e L\0,T;ll\n)) nL°°{0,T;L\n)) -. ^ G L2(0,t; (Hi(f]))')}. 

For any i9 G L^(0, r; LP(il)), we define its extension iJext by i?oxt = on [0, r] and i?ext = on (r, T]. 
It is clear that i?cxt G L''(0, T; LJ'(f2)) and the mapping ■& n> iJext is a contraction from L''(0, r; L^(i7)) 
into L«(0,r;LP(rj)). For any i) e L9(0, r; LP(17)), we define (j>r0) as the restriction on [0,t] of 0(i9cxt)- 
We infer immediately from Proposition 16.11 that ipr is continuous from L'^(0, r; L^(f2)) to L'(0, r; LP(f2)). 
Furthermore, for any 'd £ L'^(0, r; L^(f2)), we have 

\\(t>T{-&)h9{0.r;hP{n)) = \\(t)T0)h9{O.T:h^n)) = 1 1 ^l^^cxt (• , ) 1 1 l2 (o) ' < T" ^ || (/>(t5'oxt ) 1 1 L~ (0,T;L2 (O)) , 

and the previous estimates allow us to show that, for any > 0, there exists t G (0,T] such that 
(f>r maps the closed ball i3Lg(o,T;Lp(n))(0, -R'') into itself. Note that the image of i?L9(o.T;Lp(n)) (0, ^'') 
by (/) is a bounded subset of W and thus it is relatively compact in L^{0,T;LP{il)). It follows that 
the image of i3L9(o,T;Lp(n))(Oj -R'') by (fir is also relatively compact in L^(0, r; LJ'(f2)). Consequently, 
we may conclude that the problem (I3.6p - p.8l) possesses a local solution {u, z, i9) defined on [0, t] such 
that u G Wi'«(0,r;Vg(rj)), z G L°°(0, r; Hi(rj)) n Hi(0, t; L2(f])) n C0([0, r]; X,.p(17)) n L9(0, r; H2(f])), 
i, Az G L'?/2(0,r;LP(r2)) nL'?(0,T;LP/2(17)) and ^ eWr- 

We have to go back to the problem p.ip - p.3p . First we observe that g and C define a C^-diffeo- 
morphism from (0, oo) to (0, oo) and any solution of p.6p ~ p.8p provides a solution of p.ip - p.3p as soon 
as the enthalpy "i^ remains strictly positive. So we assume now that the initial enthalpy is strictly positive 
almost everywhere on fl, i.e., there exists d > such that 

5(6l"(x)) = i9°(a;) > z9 > a.e. x e fl. (6.2) 

Therefore it is possible to use the Stampacchia's truncation method and to prove a local existence result 
for the problem (|XT|) ~ (|331) . 

Theorem 6.2 (Local existence result) Assume that (IXTU)) . (|XTT|) . (IXTO)) . ((XTl) . a^d ((XTT)) 

hold. Then, for any initial data u'^ G YqIH), z° G Xq^p(f2) and i?*^ G L^(51) satisfying (j6.2p . t/iere exists 
T G (0,r] stic/i i/iai t/ie problem p.ip - p.3p admits a solution on [0,t]. 

Proof. Let (m, z, ?9) be a solution of p.6p - p.8p on [O, r] . We prove now that 

■(9(x,i) > a.e. {x,t) G Qr. 

To do so, we introduce some notations: let C ~ ^^J_ + ^ and let (/3 : [0, t] ^ M such that 

yt G [0,r] : (^(t) "^'^e-~■'^^'(^+^^""(■^'^)"L°°("))'^^ (6.3) 
We use the classical Stampacchia's truncation method. Let us define G G C^(M) satisfying 

(i) V(7 G M 3C^ > : \G'{a)\ < G^, 

(ii) G is strictly increasing on (0, oo), 
(ni) Vcr < : G{a) = 0. 

Let us define also H{a) G{s) ds for all ct G M, i?i =' — + ip and h{t) = H^-di) dx. Clearly, we 

have H G C^(IR;M) and H{a) > for all cr > 0. Furthermore, i9i(0) = -i?" + < almost everywhere 
on implies that ft,(0) = 0. Since ?9 G Wr and G H^(0,t;R), we infer that h is absolutely continuous 
and 

h{t) = [ G(??i)i9idx 
Jn 

= - / G(i9i)(div(K=V?9)+Le(M):e(M)+6l(atr(e(M))+D^i/2(^):i)+*(i)+Mi:i-95)dx 
Jn 

= - /" G'(i?i)K'=Vz9i:Vi?ida; - /" G(i9i) (Le(u):e(w)+6'(atr(e(u))+D^i?2(z):i)+*(i)+Mi:i-(y9) da:, 
Jn Jn 
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for almost every t G [0,t]. It follows from (|3.12p . p.l4[) and Cauchy-Schwarz's inequality that 

]Le{u):e{u)+a0tT{e{u)) > c^\e{u)\^ - 3a\9\\c{u)\ > ^\e{u)\'' - '^"^i^'' , 

and 

Mi:i + eB,H2iz):z > c^lip - \e\\B,H2{z):z\ > c"|ip - Cf ^ |0| (l+|z|)|i| > f jip - i£^lpEl(l+\z\^). 
Since G'{i9i) > and G{i9i) > almost everywhere and (|3.10bD and (I3.17dp hold, we get 

h{t)< f G{ih){\0\^{C+^^p^\z\^)+^)dx a.e. t e [0,t]. 

But 9 — C(i^), and reminding that (3i > 2, we have with (|3.18p that 



On the other hand, G{'di) vanishes whenever > (f, it follows that 

Ht)< / G(z9i)(|i^(C+^^|zp)+^)dx<0 a.e. t e [0,t]. 
Jn 

We may deduce that h{t) < h{0) = for all t e [0, r]. Then we infer that 

H{i}i) = a.e. {x,t) eflx (0,r), 

which implies that 

■di ^ + (fi <0 a.e. (x, t) e fl x (0, r). 
This concludes the proof. □ 



7 Global existence result 

We begin this section with some a priori estimates for the solutions of the problem (|3.6p - p.8p . As usual, 
the result relies on an energy balance combined with Gronwall's lemma. Then the global existence result 
is proved by using a contradiction argument together with the results obtained in the previous sections. 

Proposition 7.1 (Global energy estimate) Assume that ([XTU)) . (|XTT|) . (IXH)) . (|XTi)) . (IXT^ and 
(|3T7)) hold. Assume moreover that u" G Vg(f^), z" £ Xq,p(fi), -i?" G L2(f]) such that ((O)) /loWs. T/ien, 
i/iere exists a constant C > 0, depending only on ||u*'||jji(o); ll^'^llHi(f2); ll'!^'^llLi(n) ^^i^d the data such that 
for any solution {u,z,'d) of problem (|3.6p - p.8p defined on [0,t], t G (0,T], we /lawe 

Vf G [0,r] : ||K(-,f)||^,(^) + ||z(-,?)||^i(o) + m;r)hHn) < C. 

Proof. On the one hand, wc multiply (I3.6ap by ii and we integrate this expression over Q^:, with r G [0, t], 
to get 

/ (E(e(u)-z)+a6ll+Le(?i)):e(u)da;dt = / i-udxdt. (7.1) 
jQf JQf 

On the other hand, by using the definition of the subdifferential i95'(i;), we deduce from (j3.6bp that 

/ (Mi-E(e(u)-z)+D^i?i(z)+6ID^i?2(z)-i^Az):ida;di + / '^{z)dxdt = Q. (7.2) 
jQf jQf 

Adding and we obtain 



i / E(e(ti(.,f))-z(.,?)):(e(u(-,f))-z(-,f))da; + f||Vz(.,?)||2 + / Miiidxdt 

"'n JQf 

+ / Le('u):e(?i)da;di + / Hi{z{-,T))dx + / 6'(atr(e(u))+D2i?2(z):i) da;dt (7.3) 
"'Q? Jn JQf 



I ^'(i)dxdi = Co'^+ /" ^-Mdxdt. 
Jo.^ Jo.^ 
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where Cg'" = ^ J^E{c{u")-z°):{e{u°)-z°)dx + ^\\Vz°\\l,f^^^ + J^Hi{z°)dx. We integrate (IX^ over 



Qr, by taking into account the boundary conditions (I3.7p . we find 

i?(-,f)da;= / ■&"dx+ I Le(u):e(u)da;dt + / Mz:idxdt 



+ / e{atT{e{u))+'D^H2{z):z)dxdt + / ^'(i)dxdt. 
We add this last equality to (|7.3I) . and thanks to p.llap . we obtain 

i / E(eK.,f))-z(.,f)):(eK.,f))-z(.,f))dx + f||Vz(.,f)||2,(^^+c^Mk(^T)||2,(^^ 



+ / z9(-,f)da; < Cq'^ + / ^°dx + c"'\Q\+ j i-udxdt. 



Clearly there exists Ci > 0, depending only on c^, ly and c^^ such that 

Ci||u(-,r)||^,(j,)+Ci||z(.,?)||^,(f,)+ / ,?(.,?)dx<C^ + ||z9"||Li(o)+c^Mf^l+ / Mdxdt. (7.4) 
Since € G H^(0, T; L^(f2)), we may integrate by parts the last term of (|7.4p . we get 



^\H-J)\\uHn) + CilN(-,r)|||.(^) + / ^(.,f)dx < Cq"'^ + + c^'^M 

Jn 

+ IKIIcO([0,T];L2(n))||w''||L2(n) + 2Srll^llco([0,T];L2(n)) + |IKIIl2(0,T;L2(O)) + 5 / ll'"llL2(0)di, 

Jo 

which allows us to conclude with Gronwall's lemma since i!){x,t) > almost everywhere on Qr. □ 

Note that p.lSp enable us also to obtain a global estimate for the temperature. More precisely, under 
the assumptions of Proposition 17. ![ we have 

ll^llL=^(0,r;L/^i(O)) < (^5)*, (7.5) 

for any solution (m, z, 9) of problem p.ip - p.3p defined on [0, r] with t G (0, T]. 
We assume that /3i > 4. We define 

and using the notations of Section |6l we define 

Rf '^R^{R\R\ r||co([o.T];L2(o))), Rto - C(T^i?/ + ||^"||L.(o)), R"" = Tl Rto + 1- 

Then, the results of Section [6] allow us to infer that there exists r G (0,T] such that 0t- admits a fixed 
point in -Bl?(o.t;L2(si))(0i ^'')- Let us define 

f ^'^supjr G (0,T] : (/)^ admits a fixed point in BL9(o,r;L2(n))(0i ^'')} G (Oj'^l- 

It is clear that problem p.6p - p.8p admits a global solution if and only if f = T. This identity is 
established below by a contradiction argument. To do so, we assume that f G (0, T) and we choose e > 
such that f — e G (0,t). By definition of f, there exists t G {f—e,f] such that admits a fixed point 
■d = ?!'r('?) in -Bl?(o,t;L2(o))(0, -R''), i.e., the problem p.6p - p.8p admits a solution (u, z,i9) defined on 
[0,t]. We infer from the results of Section [7] that |l'i?||L°=(o,r;Li(a)) < C and ||0 = C(''^)llL9(o,T;L<i(n)) < 
Then, the definition of Rf^ and the results of Section [g] imply that = (j^ri'd) G L°°(0, t; L2(f])) with 

II'*?IIl°=(0,t;L2(O)) < Rto- 

Now let f G (0,r-f] and = ((i?'')«-f(i?^)9)l > 0. For any d G SL?(r,r+f;L2(n))(0, ^''), we 
define i?cxt as follows t9cxt = on [0, r], ?9cxt = i? on (t, t+t] and i9cxt on (r+r, T]. Clearly, we have 
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and the mapping d i^ext is a contraction on L''(t, t+t; L^(il)). Let 9 — ((dcxt)- By definition of we 
have^= =eon [0,r], = (0) on (t,t+?] andS' = C(0) = on (t+?,T]. Hence 6' G L«(0, T; L4(17)) 
and 

with = ' . By definition of (j), we get immediately that the restriction 

of 4>{'&cxt) on [0,t] coincide with 4>t{i^) = "i? and we define (/>?(i?) as the restriction of 0(??cxt) to [r, r+r]. 
Furthermore, with the estimates of Section [6l we have (f>{'dext) G L°°{0,T;L^{n)) and 

ll'/'(^^ext)||L-(0,T;L^(O)) < C {t"^ {R° , R' , M C0([0,T];L2(n)) + ll'? IIl2(0))- 

It follows that there exists r e (0,T— f], independent of r, such that admits a fixed point i? in 
T-_i_;^.L2(f2))(0, i?''). By construction of tpr, the restriction of (j){dcxt) to [0,T+r| is also a fixed point 
of (pr+T in SL«(o,r+r;L2(o)) (0, ^''). But wc may choose e 6 (0, f) such that T + T>f — e + r>f, which 
gives a contradiction with the definition of f. 

Hence we can conclude that f = T. Consequently, we deduce the following theorem: 

Theorem 7.2 (Global existence result) Assume that (jXTOt . ((3?TTt . ((3?T3)) . (|3TT4)) . (|3TT6)) and ((3Tf| 

hold. Assume moreover that /3i > 4, m° g Vo(ri); -z° G Xg^p(ri), G L^(ri) smc/i t/iat /loMs. 
T/ien f/ie problem p.6p - p.8p admits a global solution {u^z^d) such that u G W^''(0, T; VQ(r2)), z G 
L°°(0,T;Hi(rJ) nX,,p(rj)) nHi(0,T;L2(rj)), i,Az G L«/2(0,r;LP(l])) nL'?(0,T;LP/2(rj)) and G W. 
Moreover i) remains strictly positive and {u,z,9 — C,{^)) is a solution of problem p.ip - p.3p on [0,T]. 

Remark 7.3 Let us assume furthermore that there exists C^^ > and pi G (0, 1) such that 

Then we may obtain a global existence result for any /3i > max(3, ^^^3^) ■ Indeed, we can establish a 
global a priori estimate for the enthalpy gradient by using a technique inspired by Boccardo and Gallouet 
(WoG89f ). i.e., by choosing the test-function x — ^ — /'''^ some ^ > 0. By reproducing the same 

computations as in fRouKA Prop. 4-^]y '^^ may obtain an estimate ofVd in L''(0, r; L'''(f2)), independent 
of T, for any r G [l, provided that /?i > (^cLj^_2)(i-pi) ~ 5(1.^^^) since d — i. Then we consider a > 1 
such that 2^a < r and ^ >/i(^ — — for some real number /i G (0, 1). Using Gagliardo-Nirenberg's 
inequality, we may deduce that there exists Cgn > such that 

and reminding the estimate of "d in L°°(0, r; L^(il)) obtained at Proposition ]?. 1\ we infer that there exists 
a constant C > 0, depending only on the data, such that 

ll^l|L-(0,.;L»(a)) < C(l + l|Vl91|L.(0,r;L^(O)))- 

The three conditions 

1 < r < 2/ia < r, i > /i(i-i) + 1 - M with < < 1, 

allow us to choose /i — G (0, 1) and thus ct = jji = ^ + ^^"211^^ ^ (i' t) f^^'^e d — 3). It follows 

that, for any (3i > max(3, c^^i'lp^^ ) , we will obtain a global estimate of -d in L^"(0, r; L"(ri)) for any 

a G (|,|) o,nd of 9 = C,{i9) in L^^i"(0, r; L^i"(51)). Thus, with a such that (3ia > 4, we may obtain a 
global existence result by the same contradiction argument as in Section^ 
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